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We find that spin-polarized ground states emerge in nanofacets which are self-organized on SiC 
(0001) surfaces. Our large-scale density-functional calculations reveal that the nanofacet formed 
by bunching of single bilayer steps generates peculiar carbon dangling bond states localized at but 
extended along step edges. The flat-band characteristics of those C states cause either ferromagnetic 
or anti-ferromagnetic chains on covalent semiconductors. 
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Nanometer-scale structures are now accessible through 
cutting-edge fabrication technique [i| or a gift of self¬ 
organization phenomena Common expectation 

that such nanostructures break through limitations of 
current technology may rely on a fact that nano-scale 
shapes affect wave-functions of relevant electron states. 
However, a way of designing nano-shapes which deci¬ 
sively affect the electron states and therefore are intrigu¬ 
ing in science and useful in technology is lacking. 

One of the unexpected but promising examples is the 
spin-polarization or the magnetism in nanostructures 
of covalent semiconductors, that is occasionally called 
cP magnetism [M- Particular arrangements of un¬ 
paired electrons on hydrogen-covered Si (111) surface Hj 
and on Au-covered high-Miller index Si surface [mu 
are predicted to be ferrimagnetic and anti-ferromagnetic 
(AFM), respectively, by the density-functional calcula¬ 
tions. The spin-polarization on the Si surface is indeed 
evidenced by photoemissionfl^ [T^ and scanning tunnel¬ 
ing microscope [l6| experiments. 

However, such nanostructures are achieved only by so¬ 
phisticated fabrication techniques on semiconductor sur¬ 
faces. It would be free from such difficulty if the self¬ 
organization could produce desirable nanostructures by 
itself. We here report occurrence of magnetic chains on 
self-organized nanostructures of silicon carbide (SiC). 

SiC is a wide gap covalent semiconductor, existing as 
many polytypes such as 2H (wurtzite), 3C (zincblende), 
4H, 6H and so forth 17| . The 4H is the most stable with 
larger cohesive energy than others by 0.2 meV per SiC 
unit. Due to its high dielectric breakdown voltage and 
the high melting temperaturefl^. SiC is a strong can¬ 
didate for the next-generation devices in environment- 
friendly power electronics (isl. [l9| . This material shows 
peculiar nanostructures through self-organization: Dur¬ 
ing the epitaxial growth of SiC (0001) surface, the (112n) 


nanofacets are frequently observed [2rH22l | (Fig. [T]); the 
facet angle v? which is the angle between (0001) and 
(112n) is almost the same (magic angle, (p = 12° - 16°) 
and the height of the nanofacet along the (0001) direction 
is a integer or a half-integer of the unit-cell heisfht [2ll - E^ . 
The mechanism of such nanofacet formation is clarified 
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FIG. 1. (color online) Side views of the (112n) (n = 12) 
nanofacet with the facet angle (p = 15.3° on the 4H-SiC(0001) 
vicinal surface, which we have found is most stable. Green, 
yellow and purple balls indicate the Si, G and H atoms, re¬ 
spectively. The Tp‘ and Tf‘ indicate the protruded and fiat Si 
sites at the long terrace, respectively. The T®' and T^' denote 
the Si sites on the short terraces where the stacking are cubic 
and hexagonal, respectively. The Sf‘ or S®' indicates the Si 
site at the step edge on the corresponding terrace. The S^ 
and S(( denote the corresponding edge G sites. The 0i, 62 , 
and 63 indicate the bond angles at the step edge C atom. Ex¬ 
panded view around the step edge and top view of the longest 
terrace are shown in Supplementary Fig. 1 [^ . 


by our recent density-functional calculations [26|: Mi¬ 


croscopically, the nanofacet is the bunched single bilayer 
atomic steps; the balance between the step-step repul¬ 
sive energy and the surface-energy gain on the terrace is 
the reason for the formation of the particular nanofacet 
structure, strongly indicative of that such nanaofacet is 
ubiquitous on semiconductor surfaces. 


In this Letter, we show that peculiar electron states 
which are localized near but extended along step edges 
of the nanofacets of SiC exhibit either ferromagnetic 
(FM) or AFM spin polarization. Based on the density- 
functional calculations, we find that hydrogen passivation 
of the surface spontaneously induce such spin polariza¬ 
tion. 


Most calculations have been performed by the lo¬ 
cal density approximation (LDA) [ttI in the den¬ 
sity functional theory (DFT) |29|| . Key results are 
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FIG. 2. (color online) Energy bands along the FX parallel to the step edge of the stable (112n) nanofacet on 4H-SiC (0001) 
surface obtained by our LDA calculations. The Fermi level is set to be zero. The SEi and SE 2 , and also the SE 3 and SE 4 
are almost degenerate. The Kohn-Sham orbitals of the SEi, SE 2 , SE 3 and SE 4 states at the F point are shown in (b), (c), 
(d), and (e), respectively, as isovalue surfaces. The red and blue colors depict positive and negative values, respectively, of the 
Kohn-Sham orbitals. Expanded views of the Kohn-Sham orbital around the step edge are shown in Fig. [S] 


crosschecked by the generalized gradient approximation 
(GGA) [s^ to the exchange-correlation energy func¬ 
tional. We use our own real-space code named RS- 
DFT [32 h 34I c ombined with the norm-conserving pseudo¬ 
potentials [35|. We have examined effects of core correc¬ 
tions on the magnetic states and found that they 
are unimportant in the system of G and Si. The real- 
space scheme is most suitable to current multi-core mas¬ 
sively parallel architecture and allows us to perform ex¬ 
tensive calculations for thousands-atom systems. Details 
are given in Supplementary Material 


We have performed extensive total energy calculations 
for the (112n) nanofacet structures of 4H-SiG and 
obtained the most stable structure with the facet an- 
gle being 15.3° [i^, in accord with the experiments [El- 
l25l |. The structure is a nano-aggregate of bunched sin¬ 
gle bilayer steps, as shown in Fig. [T] and Supplementary 
Figure 1 (SFig. 1) [^. The bilayer atomic stacking is 
ABCBABCB • • • from the topmost surface. 


Let us start with the calculated electronic structure of 
this nanofacet (Fig. [2|). Near the Fermi level (Ep), we 
observe many states which have characters of dangling 
bonds (DBs) of surface Si atoms on the long terrace. The 
energy bands more than 1.5 eV below Ep is mainly as¬ 
cribed to the valence bands of the bulk SiC, whereas the 
bands at about 1 eV above Ep are from the bulk con¬ 
duction bands. Apart from those states, we have found 
four distinct states in the energy gap, labeled as SEi, 
SE 2 , SE 3 and SE 4 in Eig.[5])a). From the analyses of the 
Kohn-Sham (KS) orbitals of these states [Figs. [2] (b)-(e) 
and Fig. [3], we have found that they originate from the 
step-edge states: The KS orbitals of the SEi and SE 2 
states are localized near the step edges where the upper 
terraces have cubic symmetry (i.e., the stacking sequence 
from the top is either ABG or CBA; cub-terrace edge 


TABLE I. Bond angles of the step-edge C atoms in the (112n) 
nanofacet with the upper terraces being of cubic (cub) and 
hexagonal (hex) symmetry. The 61 , 62 , and 63 indicate the 
bond angles shown in Fig. [T] 


terrace 

cub 

cub 

hex 

hex 

n 

98.9 

99.3 

96.1 

96.3 

02 n 

126.8 

127.5 

129.8 

130.5 

93 n 

101.4 

101.7 

108.7 

108.6 


hereafter), whereas the orbitals of SE 3 and SE 4 states 
are near the step edges where the upper terraces have 
the hexagonal symmetry (i.e. the stacking sequence is 
either BGB or BAB; hex-terrace edge hereafter). They 
are localized near but extended along the step edges, ex¬ 
hibiting dispersion-free flat bands. 

Figure [3] clearly shows that those step-edge states have 
the character of the dangling bonds (DBs) of the three¬ 
fold coordinated edge G atoms. They are also mixed with 
the back-bond character of neighboring edge Si atoms. 
The states which have the character of DBs of the edge 
Si atoms are located around Ep, being substantially reso¬ 
nant with the terrace DB states. Larger electron affinity 
of C than Si renders the DB-character state of C to be 
located at lower energy positions being distinct from the 
Si DB-character states. 

The splitting between the SEi and SE 2 group and the 
SE 3 and SE 4 group comes from structural characteris¬ 
tics of the step edges. Table |T] shows the bond angles, 
9i, 6 * 2 , and 63 , at the step-edge G atom. We have found 
small but sizable difference in the bond angles between 
cub-terrace edge and hex-terrace edge G atoms. Apply¬ 
ing the TT-orbital axis vector analysis (POAV) [s^ to the 
three-fold coordinated G atoms where the three bonds 
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are formed by sp^ hybridization and the remaining DB 
has a character of s'^p orbital, we can deduce the amount 
of hybridization m and n from the bond angles, 9i, 62 , 
and 03 37|. The obtained values for the amounts of s 
characters in the DBs are m = 0.25 and m = 0.23 for the 
cub-terrace edge C and m = 0.13 and m = 0.13 for the 
hex-terrace edge C. The larger mixing with s-character 
renders the SEi and SE 2 states lower than SE 3 and SE 4 
states by a half eV. 

Figure[2]unequivocally show the existence of flat bands 
on semiconductor nanofacets. The flat band may gener¬ 


ally causes ferromagnetism 


and its manifesta- 
. We now find that 


tion in graphene is reported |4lL l42 
hydrogen passivation realizes the possibility of flat-band 
magnetism on the SiC nanofacet. It is of note that the 
H-passivation on semiconductor surfaces is an extremely 
common and useful fabrication technique (43l - l45j| . Table 
ITTl shows the calculated H-adsorption energies Ead 0 
for various atomic sites on the nano facet. It is found 
that Ead depends on the adsorption site. First, the H 
adsorption with the Si atoms is favorable energetically 
compared with the C atoms. As stated above, the orbital 
energy of the Si dangling bond (DB) is higher due to its 
smaller electron affinity. Before H passivation, some of 
Si DBs are doubly occupied and others are empty. The 
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FIG. 4. (color online) Energy bands [(a)-(c)] and the spin 
density [(d)-(e)] of the 4H-SiC nanofacet with terrace Si and 
some of edge Si atoms passivated by H (see text). Energy 
bands along the step edge (FX) for (a) non-magnetic, (b) 
ferromagnetic (FM) and (c) anti-ferromagnetic (AFM) states 
are shown. The Brillouin zone is a half of that in FigO The 
energy bands originated from localized states at the C atoms 
on the cnb-terrace and hex-terrace edges are marked by thick 
lines. E-p is set to be zero. The red and blue lines in (b) 
correspond to majority and minority spin states, respectively. 
The spin densities in (d) FM and (e) AFM states are shown. 
The red and blue colors in (d) and (e) represent the positive 
and negative vales for the difference between the up and down 
spin densities. Expanded views of spin densities are shown in 
SFig 2 [ 3 ^ . The color code depicting each atom is same as in 

Fig.m 


H passivation of Si DBs generates doubly occupied Si-H 
bonding states, thus leading to the larger H adsorption 
energy. Second, the H adsorption with the edge atoms is 
less favorable compared with the terrace atoms. This is 
due to the relaxation around the step edges: As repre¬ 
sented in Table HI the edge DBs have more s component 
in their characters, leading to lower orbital energies and 
less energy gain upon the H passivation. The adsorption 
energy of the Si atom on the hex-terrace edge (S^g,^) is 
relatively large. This is due to the attractive electrostatic 
interaction between Si and C in the 4H structure . 

Calculated difference in H-adsorption energies in Ta- 


FIG. 3. (color online) Expanded views of the Kohn-Sham 
orbital of the SE 3 state [Fig. [2l[d)] near the step edge. The 
top and side views are displayed in (a) and (b), respectively. 
The dashed line indicates the step edge position. The color 
code is the same as in Fig. (2] 


TABLE II. Adsorption energies of hydrogen atom Ead at var¬ 
ious Si and G sites depicted in Fig. [T]of the 4H-SiC nanofacet. 


site 

Tf 

rpbl 

-*-P 

rpbl 

cub 

rpi)l 

hex 

OSI 

'^hex 

“^hex 

qSi 

•^cub 

•^cub 

Ead 

1.71 

1.68 

1.67 

1.65 

1.56 

1.34 

1.33 

1.25 
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FIG. 5. (color online) Spin densities of the 4H-SiC nanofacet 
with terrace Si atoms and edge Si atoms passivated (a), and 
further hex-terrace edge C atoms additionally passivated (b). 
The color code is the same as in Fig. [H Expanded views are 
shown in SFig. 3 [3^ . 


bleEl indicates that the terrace Si atoms are easily ter¬ 
minated by the H adsorption, whereas the edge atoms 
are likely to survive with their DBs. The Si atom at the 
hex-terrace edge is an exception with its large adsorption 
energy explained above. Hence, when we set an energy 
criterion of 0.2 eV, it is plausible that the surface atoms 
are passivated by H except for the edge C atoms and 
the Si atoms at the cub-terrace edge. This situation is 
possible by control of temperature. 

Figure |4] shows energy bands and spin densities of 
such H-passivated nanofacet. We have found that E-p 
shifts downward and crosses the flat bands [i^. Then 
we have obtained three solutions of KS equation: The 
non-magnetic (NM), FM and AFM states. The total en¬ 
ergy of the AFM state is the lowest, being followed by the 
FM state with the energy increase of 19 meV (7.9 meV) 
per spin and then the NM state with further increase of 
83 meV (225 meV) in our LDA (GGA) calculations. In 
the NM state, Ep crosses the energy bands with the DB 
character of the hex-terrace edge C atoms (states SE 3 
and SE 4 ). The character of Si DBs at the edge is ab¬ 
sent in these bands due to the electron affinity difference 
so that the dispersion along the edge becomes smaller 
than expected from the atomic arrangement. Hence the 
exchange plus correlation splitting in the KS orbitals of 
the C DBs is 0.78 eV and 0.72 eV in the FM and AFM 
states, respectively, leading to the spin-polarized state 
in the SiG nanafacet. The AFM state in preference to 
the FM state which we have found is provisionally as¬ 
cribed to the prominent flat-band nature in the AFM 
state induced by the localization of the edge carbon DB 
(GDB) state. However, it is plausible that the AFM and 
FM states coexist in suitable temperature range. It is of 
note that the spin density [Figs. [4] (d) and (e)] is solely 
distributed along the hex-terrace edges, reflecting atom- 
scale reconstruction near edges explained above. 


The situation we have explained is characterized as the 
spin-polarization in the three-quarter filling flat bands. 
The filling by electrons is controlled by the H-passivation. 
When all the edge atoms are free from the H-passivation, 
which may be possible by increasing temperature, we 
have a situation that all GDB states are occupied and Si 
DB states are unoccupied, leading to the semiconducting 
character (not shown here). When all the Si edge atoms 
as well as the terrace Si atoms are passivated, the CDB 
states are half-occupied. We have found the AFM and 
FM states in this case in the range of the total energy of 
9.9 meV. The NM state is higher in energy by 100 meV. 
The spin density in the AFM state is shown in Fig. [5] (a). 
Interestingly, the spin density is now distributed to all 
the edges including the cub-terrace edges. This is caused 
by the filling control of the CDB states by the H passiva¬ 
tion. When the C atoms at the cub-terrace edge, which 
produce the least adsorption energy, is free from the H 
passivation, we encounter a situation in which a CDB 
state at the cube-terrace edge is half filled. We have 
found the FM and AFM states in the range of 0.1 meV 
per spin, being lower than the NM state by 141 meV. In 
this case, the spin density is distributed solely along the 
cub-terrace edges[Fig. jSKb)]. 

In summary, we have performed large-scale density- 
functional calculations on the stable 4H-SiC nanofacet 
and found that the spin polarization along the step edges 
in the nanofacet emerges. This comes from the localiza¬ 
tion of the electron states caused by the electron affinity 
difference between C and Si and by the peculiar atomic 
structures of SiC surface. We have found that the hy¬ 
drogen passivation is a powerful tool to control the spin 
polarization in the system with ubiquitous elements of 
H, C and Si. 

This work was supported by the Grants-in-Aid for sci¬ 
entific research under Contract No. 22104005 and also by 
“Computational Materials Science Initiative”, both con¬ 
ducted by MEXT, Japan. Computations were performed 
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tational Science, RIKEN and at Supercomputer Center 
in ISSP, University of Tokyo. 
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Details of Calculations 

To perform accurate large-scale total-energy electronic-structure calculations efficiently on the current and 
future architectures, we have adopted a real-space (RS) scheme[Tl [5] in the density functional theory (DFT) 
combined with the first-principles pseudo-potential scheme, and have developed a highly efficient computation 
code named RSDFT[3]. In the real-space scheme, discrete grid points are introduced in real space, and the 
Hamiltonian in the Kohn-Sham equation in the DFT is expressed as a matrix represented in terms of the 
grid-point space. Differential operators for kinetic energy are replaced by finite-difference operators with 
sufficiently high orders. Hence, the Hamiltonian matrix becomes sparse and Fast-Fourier Transformation 
(FFT) is rarely required for the Hamiltonian operation. FFT generally causes heavy communication burden 
among all the compute nodes in the conventional plane-wave basis-set scheme. As a result, the communication 
in the RSDFT code arises from the finite-difference operations, the non-local pseudo-potential operations, 
and the non-local exchange-correlation functional operations. Hence, the communication becomes local only 
among neighboring compute nodes so that the total communication cost is tremendously reduced. In fact, 
the RSDFT shows good scalability even with 80,000 compute nodes at K computer at Kobe, JapanjUl^. 

To check the validity of the results, we construct a periodic single-bilayer-step (SBS) structure without 
a long terrace. In the present calculations, the grid spacing is taken to be 0.22 A corresponding to a cutoff 
energy of 57 Ry in the plane-wave-basis-set calculations. The reliability for this grid spacing is ensured by 
the calculations with the grid spacing 0.17 A, corresponding to a cutoff energy of 100 Ry in the plane-wave- 
basis-set calculations, in the SBS structure. 

For the nanofacet formed on 4H-SiC(0001) vicinal surfaces, we treat the Si face of the (0001) surface on 
which the nanofacets are observed O [7]. To simulate the surface, we adopt the repeating slab model. We use 
an inclined slab representing the (112A^) {N = 32) vicinal surface (vicinal angle: 5.9°) and containing single¬ 
bilayer steps. The bottom surface of the slab is passivated by H atoms. Each slab consists of 7 SiC bilayers 
and separated from each other along the {112N) (TV = 32) direction by the vacuum region of 10 A. The unit 
length along the step-edge direction uq for the 1x1 lateral periodicity is taken to be 5.27 A for LDA and 
5.35 A for GGA, which is the lattice constant obtained by the total-energy calculation for the crystalline 4H- 
SiC. Brillouin-zone (BZ) integration for this slab system is performed by the four k point sampling along the 
step-edge direction: the total energy of the SBS structure varies within only 3 meV/cell when ten k points are 
used. The geometry optimization is carried out until the remaining forces become less than 50 meV /A. The 
positions of the bottom SiC bilayer passivated by H atoms are fixed throughout the geometry optimization. 
For spin-polarized calculations, in order to explore a possibility of anti-ferromagnetic couplings along step- 
edge directions, we take a double periodic cell for the step-edge direction and use two k points along the 
step-edge direction for BZ integration. We confirmed that two k points are enough to describe the magnetic 
stability by comparing the result with that obtained using eight k points. The geometry optimization in 
spin-polarized systems is performed by the same condition as the non-spin-polarized calculation. 
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Supplementary Figures 





Supplementary Figure 1: (color online) (a) Side views of the (112n) (n = 12) nanofacet with the facet angle 
If = 15.3° on the 4H-SiC(0001) vicinal surface, which we have found is most stable. Green, yellow and purple 
balls indicate the Si, C and H atoms, respectively. Expanded view around the step edge [dashed square in 
(a)] and the top view of the longest terrace are displayed in (b) and (c), respectively. 
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Supplementary Figure 2: (color online) Expanded views around the step edge for spin densities of the 4H-SiC 
nanofacet with the H-passivation except for the edge C atoms and the Si atoms at the cub-terrace edge. The 
ferromagnetic state [(a) and (b)] and anti-ferromagnetic state [(c) and (d)] are shown. The dashed lines 
indicate the step edge positions. The red and blue colors represent the positive and negative vales for the 
difference between the up and down spin densities. 
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Supplementary Figure 3: (color online) Expanded views near the step edge for spin densities of the 4H-SiC 
nanofacet with terrace Si atoms and edge Si atoms passivated [(a) and (b)] and further hex-terrace edge 
C atoms additionally passivated [(c) and (d)]. The dashed lines indicate the step edge positions. The red 
and blue colors represent the positive and negative vales for the difference between the up and down spin 
densities. 
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